Free-energy distribution of the directed polymer at high temperature 
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We study the directed polymer of length t in a random potential with fixed endpoints in dimension 
1 + 1 in the continuum and on the square lattice, by analytical and numerical methods. The universal 
regime of high temperature T is described, upon scaling 'time' t ~ T 5 j k and space x = T 3 /k (with 
k = T for the discrete model) by a continuum model with 5-function disorder correlation. Using 
the Bethe Ansatz solution for the attractive boson problem, we obtain all positive integer moments 
of the partition function. The lowest cumulants of the free energy are predicted at small time 
and found in agreement with numerics. We then obtain the exact expression at any time for the 
generating function of the free energy distribution, in terms of a Fredholm determinant. At large 
time we find that it crosses over to the Tracy Widom distribution (TW) which describes the fixed 
T infinite t limit. The exact free energy distribution is obtained for any time and compared with 
very recent results on growth and exclusion models. 



PACS numbers: 

The directed polymer (DP) in a random potential pro- 
vides the simplest example of a glass phase induced by 
quenched disorder [T] and has numerous applications, 
e.g. vortex lines [2~], domain walls [3J, biophysics g]. 
It is closely related to much studied growth models in 
the KPZ class [5], such as asymmetric exclusion pro- 
cesses (ASEP) [6J [7], and to Burgers turbulence [8]. 
It belongs to the broader class of disordered elastic 
manifolds, known to exhibit statistically scale invariant 
ground states. Within the functional RG (FRG) [5] these 
were described in a dimensional expansion by T — fixed 
points, where the ratio temperature/disorder is irrelevant 
and scales with internal size with exponent —0. 

Exact results were obtained in dimension d = 1 + 1 
PQ . Johansson proved [TOl E] that (i) the minimal en- 
ergy path of length t on a square lattice with fixed end- 
points has transverse roughness x ~ with C = 3 (ii) 
the fluctuation of the ground state energy grows as t e 
with 9 = I and its scaled distribution coincides with the 
one of the smallest eigenvalue of a hermitian random ma- 
trix, the GUE Tracy Widom (TW) distribution 12J. The 
TW distribution was found in many other related mod- 
els, polynuclear growth [15] . TASEP [6], random subse- 
quences [2][TS] and others [TMT5] , The unifying concept 
of determinantal space-time process and edge scaling was 
studied to account for such universality [T5]. An exact 
result for the space-time scaling function of the two-point 
correlator of the height in KPZ was obtained [20 . 

On the other hand, in d = 1 + 1 the model can be 
mapped onto the quantum mechanics of n attractive 
bosons in the limit n — 0, where t plays the role of (imag- 
inary) time. It can be solved with the Bethe Ansatz (BA) 
for (^-function interactions. Until now only the ground 
state energy Eo(n) was studied, i.e. the limit t — > 00 
first. Pioneering attempts at its direct analytical contin- 
uation at n = for a system of transverse size L = 00 



led to scaling behavior [2TJ [22], but not to free energy 
distribution. The possible dominance of rare events is 
also a problem of the infinite system. From continua- 
tion at fixed L, Brunet and Derrida obtained [23J the 
large deviation function for the fluctuations of the free 
energy SF ~ L 1 / 2 of the DP on the cylinder. This, how- 
ever, is different from the distribution of free energy at 
fixed t, which requires a summation over excited states. 
Also, it has been a longstanding question whether the 5- 
function model captures the low-T physics. For example, 
the FRG suggests otherwise, i.e. that some structure of 
the disorder correlator matters. In fact, Brunet's BA re- 
sult [23] for the diffusion coefficient around the cylinder, 
D ~ (kT)^ 1 / 2 , does not reproduce the expected finite 
T = limit. This remains to be reconciled with the 
standard argument of a single glass phase controlled by 
the T — fixed point, which suggests a single univer- 
sality. These issues are also outstanding for Burgers and 
KPZ growth, and related ASEP models also solvable via 
Bethe Ansatz [7]. 

In this Letter we perform the sum over excited states 
and obtain the exact expression for the free energy distri- 
bution from the Bethe Ansatz. The lowest cumulants of 
the free energy are computed at small time and checked 
with numerics. The generating function of the free en- 
ergy distribution is obtained at any time as a Fredholm 
determinant. At large time it shows that the free energy 
distribution crosses over to the Tracy Widom distribu- 
tion. The probability distribution of the free energy is 
also obtained at any time. Our study, started indepen- 
dently, parallels a recent work by Dotsenko and Klumov 
[2"5] . Although we agree on the starting sum over states 
in [25], our analysis allows to recover TW. Finally, we 
discuss the behaviour of the amplitudes as a function of 
the temperature. 

Let us recall the various definitions of the DP model. 
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(i) continuum model: the partition sum with fixed end- 
points is defined by the path integral Z = Z(0, 0, t) with: 

Z(x,y,t) = Dxe-*^Mm) 2 +v(-(r),r)] (1) 

Jx(0)=x 

for a given realization of the centered gaussian ran- 
dom potential V(x,t) of correlator V(x, t)V(x', t) ~ 
S(t — t')R(x — x'). Upon replication of ([!]), disorder av- 
eraging and Feynman-Kac formula, one finds that the 
disorder averages Z n := Z(xx,yi,t)..Z(x n ,y n ,t) satisfy 
dtZ n — —H^ ep Z n with Hamiltonian: 



i=l ij 



(2) 



and attractive interaction —R(x)/T 2 . It is known from 
FRG that to describe low-T physics one must retain some 
features of R(x), i.e. that it is a decaying function on the 
correlation scale 77. At high T, however, if one defines 



t = 2T b K~H 



(3) 



in coordinates x and i one has Z = J Dxe s , with 
S = JdtiKdix) 2 + W(x,t)} where W(x,i)W(x',i') = 
R{x-x')6(t-t') and R(x) = 2T 3 K - 1 R(T 3 k' 1 x). When 
T — > 00 one has R(z) — > 2c8(z) with c = j duR(u). 
Hence in that limit the general model ([2]), expressed in 
the coordinates x, t becomes the Lieb-Liniger (LL) model 
[ID, i.e. tHIfP = IHllL ^. of Hamiltonian: 



H. 



LL 



^ dx 2 



2c S(xi 

l<i<j<n 



(4) 



where c = — c is the interaction parameter. The LL model 
is thus the simultaneous limit T, x, t — ¥ 00 of the DP prob- 
lem with x and t fixed. It should also describe the re- 
gion where T 3 (ck) _1 ^> rf, i.e. T ^> Td ep , the crossover 
"thermal depinning" temperature, well known in vortex 
physics OEBEE]. For T > T dep the thermal fluctuations 
average out partially the disorder and R can be replaced 
by a S correlator, while for T < T dep the finite range of 
R(x) is essential for the physics of pinning. One out- 
standing question is whether, for a fixed T ^> Td ep the 
LL model describes the system all the way as it flows to 
the T = fixed point, or whether fixed but large t is a 
distinct limit from T = 0. 
(ii) discrete model: 

For the numerics we define the partition sum Z^j — 



En 



)st over all paths 7 directed along the 



diagonal on a square lattice, with only (1,0) or (0,1) 
moves, starting in (0,0) and ending in where the 

V r . s are i.i.d. random site variables. Introducing "time" 
t = i + j and space x 



2 ■ 



Z x,t 



Zi j satisfies: 



J x+i 



(5) 



with Z±fi = Sfifi. We are interested in the free energy 
F = -T\nZ with Z = Z i=o i of paths of length i re- 
turning to the origin. At T = and for a geometric 
distribution Johansson proved [TT] that the ground state 
energy Ft=o ~ e i + auji 1 / 3 with Prob{uj > — s) = F 2 (s) 
the TW distribution [JJ]. 

In the high T limit this model maps onto the con- 
tinuum one ([T]) with k = 4T and ^-function correlation 
when expressed in the variables ([3]), i.e. x = 4x/T 2 and 
I = 2i/T A . Following one checks that Z(x,t)/Z with 
Z = 2*e5' 9 * is given by the LL model Q and wc find 
that unit gaussian on the lattice corresponds to 5 = 1. 

We now use the Bethe Ansatz solution of the LL model 
Q. The moments of ([!]), expressed in x, i coordinates 
(we drop the tilde below except when stated otherwise) 
can be expressed as a quantum mechanical expectation: 



x \e- tH ^\x . 



Z n = (x 

\(x ...x \ij,)\ 2 



■ x ) , 



E 



-tE u 



HI 5 



(6) 
(7) 



i.e. all n replica start and end at xq — 0, and we used the 
resolution of the identity in terms of the eigenstates \n) 
of Hll of energies E„. The crucial observation, which 
makes the calculation tractable, is that only symmet- 
ric (i.e. bosonic) eigenstates contribute to this average. 
The eigenfunctions are superpositions of plane waves |26j 



iXp.xt 



over all permutations P 



ic sgn(xe-x k )) 



arc 



of the rapidities Xj and F[X] — Y[ n >e>e'>i ~TT 

— — v (A 

The coefficients A P = Un^eyk^ii 1 , , 

functions of the two-particle scattering phase shifts ob- 
tained from Q, and periodicity of the wavefunction re- 
quires the set of rapidities {A} to be solution to the Bethe 
equations. Major simplifications occur in this compli- 
cated equations in the attractive case c > for L = 00. 
They have complex, i.e. bound states solutions [30]. A 
general eigenstate is built by partitioning the n particles 
into a set of n s bound-states formed by rrij > 1 particles 
with n = Y)jLi n%j. The rapidities associated to a bound- 
state form a regular pattern in the complex plane which 
is called string X^ a = kj + f (j + 1 - 2a) + i5^ a . Here, 
o = l, ...,nij labels the rapidities within the string. 5 J,a 
are deviations which fall off exponentially with system 
size L. Perfect strings (i.e. with 5 = 0) are exact eigen- 
states in the limit L — > 00 for arbitrary n. Such eigen- 
states have definite momentum 

-2 c= ./-„ 2 



l m - Y.'j : '"A and 

energy E M = X)j=i( m i^ ~ fi m i( m j ~ !))• Tne ground- 
state corresponds to a single n-string with k-y — 0. The 
string-states are commonly believed to be a complete set, 
although a rigorous proof is still missing. 

To evaluate ([7]), we first obtain the string wavefunction 
at coinciding point (0 ■ • • Q|/z) = 4^(0, ..0) = n\F[X}. The 
computation of the norms \\/i\\ of string states is more 
involved, but was solved in the context of algebraic Bethe 
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Ansatz [21] (see also [35]). It reads: 
„ n2 n\{LcT° F[X] 2 



(c) n $[fc,m] 



n 



m. : 



(8) 



3[fc,m] = | I 



(ki - k.j) 2 + {m l - mj) 2 c 2 /4 
(fcj - kj) 2 + (rn, + mj) 2 c 2 /4 



l<j<7<n s 

Expressing the sum over states in ^ as all partitioning of 
n particles into n s strings and using that for L — > oo the 
string momenta rrijkj correspond to free particles [3"T] . 

Efe, -> m 3 L Jl£ we obtain & 



l.C 



^ n s !(27rc)™= 

n s — 1 v ' 



(mi , ...TTlng ) r . 



(9) 



./ = ! 



where (mi , . 



3=1 

stands for all the partitioning of 



such that Efci ' 



n with rrij > 1. 



Z = ce~T^"Z, a trivial shift in the free energy (we drop 
the hat below). This equation agrees with the one in 
[25j . although our derivation was made simpler by using 
results from algebraic Bethe ansatz. 

This formula first leads to prediction at small time. As 
in [35] we define the dimensionless parameter: 

A = {cH/l) 1 ' 3 

and z = Z/Z. Tedious calculation then yields: 

= l + \/2^A 3/2 e 2A3 (l + crf(V2A 3/2 )) 



(10) 



(11) 



In z 



32tt 



(\nz) 2 
(In^)3 C 



9-v/3 

^A 3 / 2 + (4 + 5tt - 
32 >A 3 + .. 



3?r 

Y 

32tt 



3^ 



)A 3 + 
)A 3 - 



3^3 



(12) 



The skweness of the distribution of In z is thus 7j nz w 

16v|-27 ^ 27r )l/4 A 3/4 _ 4 l/4 at smaU time _ The skewness 

for the free energy F = —Tin Z is thus 7^ = —71 z and 
negative. Fig. [T] and Fig. [2] show that the agreement 
with numerics is excellent with no free parameter. This 
is a non trivial test that the LL model is valid here and 
that the starting formula ^ is correct. 

To study any A, we avoid the explicit n = limit by 
introducing the generating function of the distribution 
P(f) of the scaled free energy F = TXf: 

g(x) = 1 + V , ; Z* = cxp(-e^-f)) (13) 

n=l 

from which P(f) is immediately extracted at A — > 00: 



lim g(x) = 9(f -x) = Prob(f > x) (14) 




FIG. 1: z 2 - 1 (4 10 6 samples) for i = 128 (triangle), i = 256 
(circle) function of i compared to formula (111 with c = 1. 



We defined — 




10 100 



FIG. 2: From top to bottom the cumulants (4 10 samples) 
(lnz) 2 (dashed line, triangle), —(In z) (solid line, circle), and 
(In a) 3 (dotted line, square) for t = 256 as compared with 
the the analytical formula ( 12 I with c = 1. 




FIG. 3: (ln 2 ) 2 7(2 4/3 A 2 ) plotted as a function of T, for in- 
creasing polymer length t. Triangles correspond to t = 4096, 
Circles to t = 256 and the dotted line to the TW variance 
0.81319... Averages are performed over 20000 samples. 
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At finite A the probability distribution of the free energy 
can also be extracted from g(x) using a Borel transform, 
equivalently described as follows. One writes (formally) 
Z = Z Z where Z = e Xu ° is a positive random vari- 
able independent from Z, with an exponential distribu- 
tion Pq(Zq) = e~ z ° (i.e. Xuq has a unit Gumbel distribu- 
tion) such that = n\Z n . The distribution P(Z) of the 
variable Z = e Xu is obtained from the cut in the grand 
canonical partition function Z(z) = Y^=o z ~ n Z n — 
as [55] : 

zP(z) = -lmZ(z + ie) = -Img(e Xx — ) (15) 

7T 7T z — ie 



with e = 0+. 

The constraint X)"=i m i = w in (|9h can then be relaxed, 
and rescaling kj — > kf/t 1 ^ 2 , it leads to: 



g(x) = l+ — f z ( n s,x) (16) 



as an expansion in the number of strings with [35] : 



7{ ,, V ff /• dkj yr (kj - kj) 2 + (m, - mj) 2 A 3 i y» ? _ 2 +Aa!m . 

l ^ , (47TA3/2)", 11 7 m , 11 (fc ! -^)2 + (m. i +m,)2A3ll e li ° 

mi,...m Tls — 1 v y j — 1 J l<z<j<n s J J J — 1 



The difficulty is the prefactor which introduces "interac- 
tions" between the strings. Let us study it in two stages: 
(i) independent string approximation. Z(l,x) can be 
computed exactly integrating over momentum: 



Z(l,x) 



dv v x l 2 

; >0 2^ 



dyAi(y)J2(-l) m e 



Xmy—vm+Xxm 



m— 1 



where, as in [25], we used that for %t[w] > 0: 



dyAi(y)e yw = e w /3 



(18) 



Rescaling v — > Aw, shifting y — >• y 
f dv v 1 ' 2 

Jv>0 



Z(l,x) 



2n 



dyAi{y + v — x) 



x we obtain: 

e Xy 



1 



after performing the sum. At large A the integration is 
only over y > and one obtains: 



lim Z(l, x) = - [ —w 3/2 Ai(w - x) 

i->oo J w>0 37T 



(19) 



Now we note that replacing in (17 1 Z(n s ,x) —> Z(l,x) n 



provides an approximation to the exact g(x) identical to 
setting the prefactor in (J7[ ) to unity: 

(20) 



■ ld {x) = exp(Z(l,x)) 



with ([19]) for A = oo and Prob md (f > x) = g m d{x). 
One easily checks that this distribution Pind(f) is the 
one obtained in Ref. [25]- Indeed, the algebraic ma- 
nipulations there are equivalent to setting the prefactor 
to unity at large A [37 . However this distribution has 
skewness 71 = 0.96029, incompatible with our numer- 
ics which shows instead for all A a negative skewness 7^ 



bounded by (minus) the TW skewness 71 = -0.224084... 
Although one checks that it reproduces the leading tail 
for / —> —00 of the TW distribution, it differs from it. 
As we now show, including interactions between strings 
leads to TW. 

(ii) exact result for the generating function at any time: 
We now derive an expression of g(x) valid for any A, 
in terms of a Fredholm determinant. Using the identity: 



det[- 



1 



n 



i(ki — kj)X 3 / 2 + (mj + m,j) 

2\3 



(h — kj) 2 + (mi — mj) X 
(ki - kj) 2 + (mi + mj) 2 X 3 



n 



2to,- 



(21) 
(22) 



and manipulations as above, starting from ( 17) one finds 



Z(n s ,x) 



Vi>0 ,._ 



J\dv t det[K x (vi,Vj)] (23) 



i=l 



g(x) = Det[l + P K X P Q ] 



(24) 



where Det is a Fredholm determinant (FD) defined with 
integration on the real positive axis J v>0 , i.e. here and 
below we define P s the projector on [s, +00 [. The kernel 
K x (v, v') — <$> x (v + v , v — v ), where we have defined the 
function: 



dk 

2^ 



dyAi(y 



x + u) 



pXy—ikw 



1 



,Xy 



(25) 



These formula generate the small A expansion but they 

are valid for all A. 

(iii) free energy distribution in the large time limit 
For large A one can replace 1 ^ e " B — > 0(y). Then one 
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obtains for A 



-oo: 



Prob(f > x) = g(x) = det(l + P„*KP_*) (26) 
f dk 

K(v, v') = - —dyAi(y + k 2 + v + v')e~ ik{v - v ) 

Jy>0 27r 

where all integrals in the FD are for f . ,„. One re- 

G Ji?>— x/2 

covers in particular Z(l,x) = TrK x = J v> _ x / 2 K(v,v) 



which yields ( 19 ) above. We can now use the following 
identity between Airy functions [55] : 



dkAi(k 2 +v + v , )e lk( - v - v ' ) = 2 2 ' 3 T:Ai{2 1 ' 3 v)Ai{2 1 ' 3 v') 

(27) 



which immediately implies that: 

K(v,v r ) ^^^Km^v^v 1 ) 



(28) 



where K Ai (v,v') = (Ai(v)Ai'(v') - Ai'(v)Ai(v'))/(v - 
v') = J y>() Ai(v + y)Ai(v' + y) is the Airy kernel. Upon 

rescaling v, v' by a factor 2 -1 / 3 we obtain: 

Prob(f >x = -2 2 / 3 s) = Det(l - P s K Al P s ) = F 2 {s) 

(29) 

i.e. the Tracy Widom distribution. Hence in the large 
time limit one recovers the TW distribution. 
(iv) free energy distribution for any time 
We now extract the free energy distribution for any 



time. We use ( 15 ) which expresses the distribution of 
In Z as a convolution, i.e. In Z = Xuq + Am where A«o 
is a unit Gumbel independent random variable and the 
distribution of the variable Xu = In Z is obtained as: 

p(u) = -?-(£>ei(l + P KP ) - Det(l + P a K*P )) 

2tTT 

/dk e \y-ik(v-v') 
-dyA l {y + k 2 +v + v') eXu _ K _ eXy 

(30) 

where all / u>0 and * denotes complex conjugation. Using 
l/(x— ie) = PV\ +iir5(x), the complex kernel is written 
K = K 1 + iK 2 where, using (27) one finds: 



Ki (v,v') = 25 PV / dy 



Ai(2^v + y))Ai(2iv' + y)) 



„Xu-2 2 / 3 \y _ I 



(31) 



and 



K 2 (v,v') = ^P Al (2^v + 2- 2 / 3 u,2^v> + 2~ 2 /*u) 

(32) 

where Pai(v, v') — Ai(v)Ai(v') is a rank one projector. 
The latter property implies that p(u) is a linear function 
of K 2 hence: 

p(u) = Det(l + P {Kx + X K 2 )P ) - Det(l + PqK^) 

7T 

(33) 



which is our final expression [56] . 

We can now compare our result with the very recent 
works 39J on KPZ growth with the narrow wedge initial 
condition, to which our work also applies (see also [ID]). 
The correspondence reads that A ^ z h = \nZ, 2v = T/n 
and DX 2 KPZ = c/k 2 . There the distribution of h was 
obtained, which translated here yields (up to an additive 
constant) \nZ — y£t where 7 = 2 2 / 3 A and the distri- 
bution of £ t becomes identical to TW at large t. Upon 
rescaling of v,v' by 2 -1 / 3 our result (34) can also be 
rewritten as: 



2 2 l z V {u)=Det[\-P 2 _ lu (B t 
-Det\l-P 2 B t P 2 1 

L 2 3 u 2~ 3 u' 



Pm)P. 



2 3 u' 



(34) 



hence the results 



where B t is the kernel defined in 
coincide. 

Let us close by discussing the temperature dependence 
for experimentally relevant models, e.g. either R(u) with 
a finite range correlation, or a discrete model. For any 

fixed T one expects AF = F^ 1/2 = A{T)t e at large t, 
with 9 = |. Concerning the amplitude A(T) it is clear 
that the (5-function model reproduces only its high T 
behaviour. Indeed, here we found AF = Tf(t), with 
f(i) - t 1 / 4 at small t from |Sj) and f(t) - t e from 
our large A analysis, 
tude A(T) 



Hence large A yields the ampli- 
K l/3y-2/3 anc j ^ n j g can Qn iy ^ e interpreted 

as a high T limiting behavior, i.e there is no way the 5- 
function model can predict the amplitude for the distinct 
T < Tdep regime, where it crosses over to a constant A(0). 
This is illustrated in Fig. [i]where a 2 = (lnz) 2 7(2 4/3 A 2 ) 
is plotted as a function of T, for increasing t. The fixed 
T and large t behavior is a 2 ~ A(T) 2 T^ 3 k~ 2 ^, hence at 
low T it behaves, for the discrete model as ~ ^4(0) 2 T 2 / 3 
with a non-universal prefactor. At high T we know from 
the small A prediction ( 12 ) with 5=1 that it behaves as 
~ T 2 / 3 ^ 1 / 6 . For intermediate T a plateau is thus pre- 
dicted to develop. Its approach from above is described 
by the large A limit of the (universal) crossover function 
computed here. The low T behaviour below the plateau 
is out of reach of the (5-function model. The value of the 
plateau should equal the variance of the TW distribution 
(j 2 = o\ w = 0.81319... However one sees that the con- 
vergence is slow and requires very large polymer lengths 

m 

Similarly one surmises that Ax ~ T 3 g(t/T 5 ) with 
g(y) ~ y 1 ! 2 at small y and g(y) = y^ at large y interpolat- 
ing between thermal diffusion x ~ y/Tt and x ~ B(T)t 2 / 3 
with B(T) - (kT)- 1 / 3 at high T, while for T < T depi 
B(T) « B(T = 0). A similar interpretation of Brunet's 
result for the winding D ~ (kT) -1 / 2 can be given. Note 
that if the exponents were to assume their Flory values, 
Qp = fj, 9p = then high and low T regimes would 
merge without need for a plateau (i.e. A(T) is constant 
at high T), and this is indeed the case in the mean field 
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method [H]. This is not the case here, and this is be- 
cause in the high T regime more typical paths contribute, 
AF = 0(T), while the low T problem is dominated by 
the lowest energy path, AF = 0(1). 

To conclude we have obtained the distribution of the 
free energy of directed polymers in the high T regime 
described by the attractive Lieb Liniger model, from the 
Bethe Ansatz. It becomes identical to the Tracy Widom 
distribution at large time, although the amplitudes ex- 
hibit a distinct low temperature behavior. 
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